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Abstract 

We model the physical behaviour at the surface of a relativistic radiating star in the 
strong gravity limit. The spacetime in the interior is taken to be spherically symmet¬ 
rical and shear-free. The heat conduction in the interior of the star is governed by the 
geodesic motion of fluid particles and a nonvanishing radially directed heat flux. The 
local atmosphere in the exterior region is a two-component system consisting of stan¬ 
dard pressureless (null) radiation and an additional null fluid with nonzero pressure 
and constant energy density. We analyse the generalised junction condition for the 
matter and gravitational variables on the stellar surface and generate an exact solu¬ 
tion. We investigate the effect of the exterior energy density on the temporal evolution 
of the radiating fluid pressure, luminosty, gravitational redshift and mass flow at the 
boundary of the star. The influence of the density on the rate of gravitational collapse 
is also probed and the strong, dominant and weak energy conditions are also tested. 
We show that the presence of the additional null fluid has a signihcant effect on the 
dynamical evolution of the star. 

Keywords: two-fluid atmospheres; compact radiating stars; relativistic astrophysics 

1 Introduction 

A compact star is formed when a massive star (M* > 8Mq) breaks away from a state 
of hydrostatic equilibrium and collapses under its own gravity. This situation usually 
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arises when the massive star has reached the end of the first phase of its evolution. 
The stellar object that is formed at the end of the gravitational collapse process is 
very dense and small in size with a typical radius of the order of (7 — 20)fcm for ultra¬ 
compact quark-gluon stars, strange stars and pulsars (including unmagnetized neutron 
stars). As a result of the high mass density in the resulting conhguration, the gravi¬ 
tational field throughout the interior is significantly strong and the star has effectively 
transitioned into the so called strong gravity regime. At this point the object is a 
far more enhanced self gravitating system and is termed a relativistic star since the 
theory of general relativity would be required if we are to construct a plausible and 
realistic model for its global dynamics. For detailed and comprehensive reviews on the 
theory of nonadiabatic gravitational collapse and relativistic stars the reader is referred 
to the works of Oppenheimer and Snyder-, Penrose^, Glendenning-, Strauman- and 
Shapiro and Teukolsky-, among others. During the contraction of the massive star 
and the evolution of the compact object, gravitational binding energy is converted into 
heat energy which is used for ionization and dissociation in the dense interior. The 
excess heat energy is dissipated as radiation across the stellar surface. This process is 
described adequately for a relativistic star by constructing the junction conditions at 
the surface. These equations relate the interior and exterior matter variables as well 
as the respective geometries. 

It was in the pioneering work of Oppenheimer and Snyder- that the signihcant 
role of general relativity in modeling the nonadiabatic gravitational collapse of stellar 
bodies was established. The dynamics was studied by considering the contraction of 
a spherically symmetric dust cloud and it was realised that relativistic gravity could 
adequately describe the outflow of heat energy in the form of pressureless radiation 
which is often referred to as null dust. While this seminal investigation provided an 
initial framework in which to construct models for the internal evolution of the dissi¬ 
pating object, it was not completely clear what the conditions were at the boundary or 
surface. Moreover it was also not clear as to what the appropriate geometry of the ex¬ 
terior region was, i.e. there was no known exact solution to the Einstein field equations 
that described the gravitational held of a radiating star. Interestingly, it was some 
time later that the hrst solution that describes this scenario was obtained, and this has 
since then become a crucial and essential ingredient when studying stellar like objects 
with heat how and strong gravitational helds in astrophysics. The Vaidya^^ metric is 
the most well known and widely used solution in general relativity that describes the 
exterior gravitational held of a dense object exhibiting outward radial heat how. It 
dehnes outgoing null radiation and is written in terms of the mass of the radiating 
body. This result provided a marked advancement in the modeling process and it was 
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eventually when the so called junction conditions for radiating stars was generated by 
Santos^ that general relativity provided a satisfactory framework in which these objects 
could be studied in greater detail. The Santos junction conditions are physically very 
meaningful, for localised astrophysical objects because they reveal that the pressure of 
the radiating stellar fluid is not zero at the boundary but rather that the pressure is 
proportional to the magnitude of the heat flux. This seemingly simple mathematical 
relation has to be solved as a highly nonlinear differential equation; the solutions of 
which provide the evolution of the gravitational held potentials and in turn the be¬ 
haviour of the matter held. However, it must be noted that this standard framework 
only describes the emission of pressureless null radiation (photons) into the exterior 
region of the dissipating relativistic star. It has not been used to probe the outhow of 
any other type of observable radiation or elementary particles like neutrinos, which are 
thought to be signihcant carriers of heat energy in stars and released by particle pro¬ 
duction processes at the stellar surface. Recently, Maharaj et al— have generalised the 
Santos junction condition by matching the interior geometry of a spacetime containing 
a shear-free heat conducting stellar huid to the geometry of an exterior region that is 
described by the generalised Vaidya metric which contains an additional Type II null 
fluid. This new result provides a greater array of possibilities with regards to the mod¬ 
eling of relativistic objects in astrophysics as it depicts a more general exterior region 
for a radiating star, which is comprised of a two-fluid system: a combination of the 
standard null radiation as in the case of the Santos framework, and an additional fluid 
distribution which is more general and can be taken to be another form of radiation or 
better still, a held of particles such as neutrinos, as mentioned above. An interesting 
feature of this generalised junction condition is that the radiating huid pressure at the 
boundary is not only coupled to the internal heat hux but also to the non-vanishing 
energy density of the Type II null huid. This conveys a direct relationship between the 
evolution of the interior and exterior matter and gravitational helds, at the boundary 
of the star and consequently may yield physical behaviour that is far diherent and 
perhaps more realistic from that which arises in the standard scenario. 

A substantial amount of work on relativistic radiating stars has been carried out in 
the standard Santos framework. In the investigations of de Oliveira et al— and Maharaj 
and Govender— , the Santos junction conditions were generalised to include the ehects 
of an electromagnetic held and shearing anisotropic stresses during dissipative stellar 
collapse. More recently, analytical models for shear-free nonadiabatic collapse in the 
presence of electric charge were obtained by Pinheiro and Chan—. The inhuence of 
pressure anisotropy, shear and bulk viscosity on the density, mass, luminosity and ef¬ 
fective adiabatic index of an 8Mq contracting radiating star was studied by Chan— >1^. 
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These results were later extended by Pinheiro and Chan—. Misthry et al— generated 
nonlinear exact models in the shear-free regime for relativistic stars with heat flow, 
using a group theoretic approach that involves the Lie symmetry analysis of the San¬ 
tos junction condition. Following this, Abebe et al—, employing the same technique, 
investigated the behaviour of radiating stars in conformally flat spacetime manifolds. 
Abebe et al—^ used Lie analysis to hnd radiating Euclidean stars with an equation of 
state. Collapse models with internal pressure isotropy and vanishing Weyl stresses were 
also probed by Maharaj and Govender— . They investigated the dynamical stability of 
the dissipating stellar fluid and demonstrated that the configuration was more unstable 
close to the centre than in the outer regions. It is also well understood that the thermal 
evolution of the radiating fluid is crucial in any stellar model and the precise role of the 
relaxation and mean collision time were analysed by Martinez^, Herrera and Santo&^i 
and Govender et al—. These ideas were further exploited in the more recent works 
of Naidu et al—, Naidu and Govender— and Maharaj et al—. Attempts to model 
radiating stellar matter with a more realistic form, have been made by imposing either 
a barotropic or polytropic equation of state for the fluid distribution. Wagh et al— 
implemented a linear equation of state for matter evolving in shear-free spacetimes and 
Goswami and Joshi— studied the gravitational collapse of an isentropic perfect fluid 
with a linear equation of state. 

It is noteworthy, at this point, to mention that not much work has been done in 
describing dissipating stellar bodies in the more general scenario when the exterior 
region is a two-fluid system and the Santos junction condition is generalised. Since 
the appearance of the Maharaj et al— result in the literature, no exact solution and/or 
corresponding physical model has been presented. However, it should be noted that 
the idea of a Type 11 null fluid existing in the local exterior of a radiating star has been 
explored previously in other contexts and without there being any direct connection 
between the interior and exterior. The physical properties of the generalised Vaidya 
spacetime containing a Type 11 null fluid were studied extensively by Husain—, for 
which exact solutions were obtained strictly for the exterior fluid. In a subsequent 
investigation Wang and Wu— further extended theses ideas and generated classes of 
solutions for spherically symmetric geometries. These general results contain the well 
known monopole solution, the de Sitter and anti-de Sitter solutions, the charged Vaidya 
models and the radiating dyon solution. The conditions for physically reasonable en¬ 
ergy transport mechanisms in the generalised Vaidya spacetime were explored by Krisch 
and Glass^ and Yang—. Glass and Krisch— used the notion of the null fluid to de¬ 
scribe the qualitative features of a localised string distribution in 4-dimensions, and 
generated analytical solutions for null string fluids exhibiting pressure isotropy in the 


4 


limit of diffusion. This realisation that a Type II null fluid can be used to describe 
a diffusing string on cosmological scales, in particular, supports the earlier works of 
Vilenkin—, Percival— and Percival and Strun&2i. 

This paper is organised as follows: In section 2 we present the basic theory for 
relativistic radiating stellar models with spherically symmetric shear-free interiors and 
two-fluid exteriors. The generalised junction condition is then integrated as a non¬ 
linear second order differential equation at the boundary, and an exact solution for 
constant null fluid energy density is generated. The luminosities and gravitational 
redshift are then defined for the case when the null fluid is present on the outside. 
Finally, some qualitative features of the null fluid are briefly described. A detailed 
physical analysis is carried out in section 3. The evolution of the radiating fluid and 
the gravitational variables at the surface, are probed and the standard and generalised 
models are compared in order to establish the impact of the exterior null fluid. This is 
achieved by constructing temporal prohles of the radiating fluid pressure, the surface 
and asymptotic luminosities, and the gravitational redshift using the exact solution due 
to Thirukkanesh and Maharaj— and the solution presented in this work. The strong, 
dominant and weak energy conditions, mass flow, and rate of gravitational collapse for 
the two scenarios are finally investigated. In section 4 we discuss our results and make 
concluding remarks. 

2 Constructing the model 

2.1 The interior and exterior geometry 

The dynamics of the gravitational field in the stellar interior, in the absence of shearing 
stresses, is governed by the spacetime line element 

= —A^(r, t)dt^ + t)[dr‘^ + r‘^{d6‘^ + sin^ dd0^)], (1) 

where A{r, t) and i?(r, t) are the relativistic gravitational potentials. The matter inside 
the star is defined by a relativistic fluid with heat conduction in the form of the energy 
momentum tensor 

Tab = (/^ + P)UaUb + PQab + QaUb + qbUa- (2) 

Here /i and p are the fluid energy density and isotropic pressure respectively. The 
components Qab represent the metric tensor field and u and q, are the fluid four-velocity 
and the radial heat flux, respectively. The Einstein field equations = T~^ for the 
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interior may be written as 
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where dots and primes denote differentiation with respect to coordinate time t and 
radial distance r respectively. 

In the local region outside the star the dynamics of the gravitational held may be 
described by the generalised Vaidya outgoing radiation metric which has the following 
form 

ds^ = — ^ dv"^ — 2dvdr + r^{d6‘^ + sin^ OdcjA), (4) 

where m{v, r) represents the mass how function at the surface, and is related to the 
gravitational energy within a given radius r. The characteristic feature about the metric 
()Tj) is that the mass function also contains a spatial dependence in the radial direction 
during dissipation; this is signihcantly diherent form the standard scenario in which 
the mass at the boundary only has a time depenhance. Husain— and Wang and Wu— 
have shown that an energy momentum tensor that is consistent with (jl]) is 


TA — £^Jb + (p + -P) i}a^b + h^a) + Pdabj 


( 5 ) 


which is a superposition of null radiation and an arbitrary null huid. Here e is the 
energy density of the photon radiation, and p and P are the energy density and pressure 
of the additional null fluid, respectively. In general, represents a Type H fluid 
as dehned by Hawking and Ellis^. This additional fluid component has also been 
interpreted by Husain—, Wang and Wu— , Glass and Krisch— and Krisch and Glass^ 
as a string fluid distribution in four dimensions. The null vector 1 is a double null 
eigenvector of the energy momentum tensor ([5]) and the vector n is normal to the 
spatial hypersurface. For the spacetime metric (jl]) and the energy momentum tensor 
(jSj) the Einstein field equations for the external local two-fluid stellar atmosphere can 
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be written as 
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It is interesting to note that in the standard framework, the Santos^ junction con¬ 
dition tells us that on the boundary of the star the pressure of the stellar fluid is 
proportional to the magnitude of the heat flux 

(p)s = {qB)T.- (7) 

This condition has recently been extended to incorporate the generalised Vaidya ra¬ 
diating solution. Maharaj et al— have shown from hrst principles that the junction 
condition (describing the dissipation) for the metric (jl]) can be written in the form 

(p)s = {qB - p)s. (8) 

Note here, that this result is more general and contains the Santos equation ([7]) as a 
special case when p = 0, i.e. when the additional null fluid is absent and the outside 
contains only pure radiation. It is also important to observe that the generalised 
boundary equation ([8]) is new for spherically symmetric nonadiabatic stellar models and 
consequently have no reported exact solutions. Furthermore, the new result is quite 
extensive and may be applicable to more physically realistic astrophysical scenarios 
as it includes the effect of the null fluid energy density. As mentioned earlier, the 
null fluid and its nonvanishing energy density have been explored in the context of 
cosmological and more localised string fluid distributions in four dimensions. It may 
also be appropriate for the description of neutrino outflow from compact relativistic 
stellar objects in which nonadiabatic and particle production processes prevail. In view 
of this it is important to emphasize that the generalised null fluid models should allow 
for signihcant improvement on the results obtained by Glass^, following the seminal 
treatment by Misner— . 

2.2 Qualitative features of the exterior null fluid 

For a better understanding of the impact of the additional null fluid in the exterior of 
the radiating body, it is worthwhile to establish what the general features and qualities 
of the fluid are. It is clearly evident from the exterior held equations fl^ and fl6cl) 
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that an explicit relationship between the energy density p and the pressnre P can be 


obtained. By simply differentiating throngh fl^ with respect the exterior radial coor¬ 
dinate r and using we see that 



(9) 


This equation indicates that there exists, in general, a linear connection between the 
pressure and the energy density. Moreover, (|9]) also suggests that given a form for the 
density, we can determine both the magnitude and the sign of the pressure. The sign 


of the pressure is crucial as it reveals more about the nature of the null fluid and the 


interaction of the fluid with the gravitational field. If an arbitrary functional form is 
prescribed for the density then it is more difficult to determine the sign of the pressure, 
as this would depend on the actual expressions for the terms in brackets. However, 
if as in the case of the model presented in this work, the density at the surface is 
taken to be constant then it is straightforward to obtain the signature of P since for 
physically meaningful applications the energy density has to be strictly positive (p > 0). 
With this in mind it turns out that for a constant null fluid energy density at the 
surface, the pressure P is always negative. This is consistent with the ideas exploited 
in cosmological scenarios where a fluid characterised by the cosmological constant or as 
dark energy or a scalar field, has to have a negative pressure in order to counteract the 
attractive force of gravity induced by matter and large scale structure in the universe. 
In a more localised astrophysical setting as in this study, the interpretation may be 
somewhat different. 

2.3 Luminosity and gravitational redshift 

In our new generalised framework, in which the exterior region of the radiating body 
consists of both standard null radiation (photons) and an additional Type II null fluid, 
we expect that the null fluid should have some effect on the observable and measurable 
properties of the radiating fluid at the surface. It is quite evident from the generalised 
junction condition ([8]) that the fluid pressure at the boundary now depends on the 
energy density p of the null fluid. Moreover, it has also been well established by Bon- 
nor et al— that the luminosity and redshift of a dissipating star in general relativity 
depend on the fluid pressure; hence one expects the null fluid to in some way influence 
these quantities at the surface. 


We aim now to generalise the definitions for the fluid luminosity and gravitational 
redshift in the context of our new model. The null radiation emitted by a relativistic 


object with heat flow, experiences a gravitational redshift (change in energy) due to 
its motion in the gravitational held. This redshift is generally written as 


= 



( 10 ) 


where is the luminosity of the radiating fluid at the surface and L^o 


is the luminosity 


of the radiating fluid as seen by a stationary observer at an infinite distance away from 
the object. These are given respectively, by 
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In the above m is the mass flow across the boundary E, r is the timelike coordinate 
on the boundary, and v is the timelike coordinate in the exterior region described by 
the generalised Vaidya spacetime (jl]). (For details see Maharaj et al—). The term ^ is 
given in terms of the interior gravitational potentials by 


dv 


The derivative is written as 

dr 





The mass flow is given by 
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Differentiating flT^ results in the expression 

fdm^ _fr^B^ r^BBB r^BB^A r^B'B' B'^b\ 


( 13 ) 


( 14 ) 


( 15 ) 


( 16 ) 


9 












Then using the generalised junction condition ([8]) as well as (l3b|) in flT6|) gives 
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dt 
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( 17 ) 


which is the generalised result for the mass flow rate when the null fluid is present in 
the exterior with the generalised Vaidya metric with m = m(n, r). When the null fluid 
is absent (p = 0) we then get 

^ = -1 (^pr^B^B + qr^AB‘^{B + rB')^ , (18) 


which corresponds to the mass flow rate when only null dust is present in the exterior 
with the standard Vaidya metric with m = m{v). Hence we have shown that the mass 
flow rate is fundamentally affected by the 2-fluid atmosphere with the generalised 
Vaidya metric since the null fluid density p is nonzero. 

Using the derivatives flTT)) . flTT)) . and flT^ in the dehnitions flTT]) and flT^ we generate 
the following expressions for the surface luminosity Ls and the asymptotic luminosity 
Loo, respectively 


Lt. 
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OD 
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(19a) 

(19b) 


Equations fll9al) and fll9bl) are the generalisations of the surface and asymptotic lumi¬ 
nosities when the null fluid is present in the exterior region. It is clear that the null 
fluid should have a marked impact on the values of these quantities, as the energy 
density p appears explicitly through the boundary pressure ps = (<?-B “ p)s, and the 
gravitational potentials A and B, in the above forms. Furthermore L^ and L^o differ 
by a factor of j • We observe, once again, that when the null fluid is absent 

in the exterior (p = 0), fll9ap and fll9&P reduce to the following forms. 


Ls = -r'^B'^p, 
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which are just the classical definitions for the luminosities when there is only radiation 
in the outside of the dissipating body. We emphasise the point that the luminosities 
Ls and L^o are fundamentally charged by the 2-fluid atmosphere with the generalised 
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Vaidya metric with m = m{v,r) since the null fluid energy density affects the space- 
time geometry. Investigating the spatial and temporal behaviour of and L^o is of 
crucial importance when modeling dense stars in astrophysics. For a particular model 
of a spherically symmetrical relativistic star with heat flow, Tewari^-^ used his exact 
solution at the stellar boundary to compute Ls and Loo- He demonstrated that the 
luminosity at infinity could be as much as 2.61089 x 10^® ergs/sec and that the surface 
luminosity was of the order of 9.8 x 10^® ergs/sec. Sarwe and Tikekar— analysed the 
evolution of the asymptotic luminosity and used it to determine the timescale for which 
a 3.24Mq compact star with radius in the order of 17km collapses to form a black hole. 


Utilising the system fll9p and the definition fllUp we may write down the form for 
the gravitational redshift of the radiating fluid at the surface as 


zt. 



kikV -1 
« ) 


( 21 ) 


Although it appears that the structural form of the redshift has not changed, the mag¬ 
nitude of zj] will be different since the potentials A and B are obtained by solving 
equation (|8]) (which when written in expanded form is a second order nonlinear differ¬ 
ential equation that describes the generalised two-fluid system in the exterior) at the 
boundary. Computing the gravitational redshift is useful when modeling relativistic 
objects with heat flow as the form and the magnitude of z-£ can be used to give a 
measure of compactness or denseness of the object; that is to say, measuring the red¬ 
shift determines the compactness of the relativistic body. Hladik and Stuchlik— have 
studied the gravitational redshift of photons and neutrinos radiated by neutron stars 
and quark stars in the braneworld scenario. They pointed out that when observational 
data for the photon and neutrino surface redshift is known, their model can be used 
to determine the critical parameters that describe the compactness of the star under 
investigation. 


2.4 Generating an exact solution 

We now turn our attention to the process of generating an exact solution to the new 
boundary condition (jH]). Considering the high degree of nonlinearity and other mathe¬ 
matical complexities that surround equation (|H]) it is convenient at this point to focus 
on relativistic stellar models in which the fluid particles move along null geodesics 
from the core and up through to the stellar surface where the radiation is lost to 
the exterior. Such models were investigated by Rajah and Maharaj— , Govender and 
Thirukkanesh— and Thirukkanesh and Maharaj—. In particular, the Govender and 
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Thirukkanesh— model investigated the behaviour of a radiating star by including a 
cosmological constant term in the Einstein held equations for the interior. It must be 
pointed out, however, that although their resulting junction condition has the same 
structure as in this study, their model only accounts for the presence of pure null radi¬ 
ation in the exterior. The motivation for this idea is that the dissipating star may be 
immersed in an ambient environment that contains dark energy which is not connected 
in any way to the interior of the body. We do not follow the same line of thought but 
treat the null huid distribution in the most general way, and regard it as more likely 
to be sourced from within the star during the radiating phase. Consequently the fluid 
may be a form of high energy radiation or a held of elementary particles such as neu¬ 
trinos or even electron-positron pairs as in the case of strongly magnetized neutron 
stars. The studies mentioned above, were carried out in the framework of the standard 
Santos junction condition p = qB. It is our aim now to extend these models by utilising 
the new generalised junction condition ([8]), and make use of the fact that the stellar 
atmosphere is now a well dehned two-huid system. For our investigations, we consider 
the null huid energy density p to be constant on the stellar surface. This assumption 
is not physically unreasonable as this situation corresponds to the dihusion of the null 
huid in the stellar exterior as the relativistic star approaches a state of hydrostatic 
equilibrium and enters the early stages of thermal cooling. 

For geodesic motion of huid particles we have the gravitational potential A = 1. 
The condition of pressure isotropy in the absence of shearing huid stresses admits the 
following analytical form 


Ci(f)r2 - C2(f)’ 

for the gravitational potential B. It is also important to recall that the above form 
corresponds to conformally hat tidal gravitational ehects and was hrst obtained by 
Banerjee et al—. Here d is an arbitrary constant and Ci{t) and C 2 {t) are functions of 
integration which have to be determined in order to complete the exact solution. The 
huid pressure and radial heat hux reduce to the following equations 
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from (|3b|l and (|M|l respectively. With the above system (|23|1 and the form (|22l) the 
generalised junction condition given by (E]) can be fully expanded as 


- 4db(ciC2 - ciC2)(cib^ - C 2 ) - 4ciC2(ci6^ - 02 )^ - 2d‘^{cib'^ - C2)(ci6^ - C 2 ) 
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+ 5(f{cib‘^ — C 2 Y — p(f{cib‘^ — C 2 Y = 0, (24) 

where we have taken r = r-£ = b{= constant), on the stellar surface. The value 
= & is the actual radius of the stellar distribution and for compact and ultracompact 
relativistic bodies like neutron stars and pulsars, and strange or quark stars, is of the 
order of 7 — SOfcms. It is important to note the appearance of the additional term that 
arises due to the presence of the null fluid energy density p 7 ^ 0. This parameter is 
taken to be constant on the stellar boundary and as mentioned earlier, is an appropriate 
and reasonable assumption. In the limit when the density p goes to zero, (for which 
the stellar exterior is composed only of pure radiation), we regain the corresponding 
equation obtained by Thirukkanesh and Maharaj— in the context of the standard 
framework. That equation has been studied thoroughly and one of the earliest known 
solutions obtained was presented by Kolassis et al—. Our new equation fl2T)) now has an 
added degree of complexity as a consequence of the null fluid component and analysing 
it is more involved than in previous cases. 

We realise that in order to integrate fl24)) . the following transformation can be 
utilised 

u{t) = Cib^ — C 2 - (25) 

Then equation fl2T)) may be rewritten as 


Abdu^ci + 4(-u^ — bdu)uci — Ab‘^u^c\ = \{2uu — bii^) + . (26) 


Equation fl26|) is a Riccati equation in ci but is still difficult to solve in general. If we 
let M = a (constant) then becomes 


a bod 

+ r/' - = 46 ^' 


We now make use of the transformation 


Cl = - 


d w 
b ta’ 


(27) 


(28) 


where w{t) is an arbitrary function. Then equation fE7)) becomes 


w + = 0, (29) 

bd 4 

which is a second order linear ordinary differential equation with constant coefficients. 
The general solution to equation (|29|) is given by 


w(i) ^ 9i{t) exp [Aii] + g 2 {t) exp [A 2 i] . 


(30) 


In the above solution gi{t) and g 2 {t) are functions of integration and 



( 31 ) 
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Then the functions Ci(t) and C2(t) become 


Ci(t) 


C2(t) 
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'giXi exp(Aif) ^faAa exp(A 2 t)' 

V h) 

5fi exp(Aif)5f2 exp(A 2 t) 


-{bd) 


giXi exp(Ait) + g 2\2 exp(A 2 t) 

5(iexp(Ait) + 5(2exp(A2t) 


— a. 


Consequently the gravitational potential B has the form 




-db 
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glAi exp(Aif)+g2A2 exp(A2p 
exp(Aip+g2 exp(A2t) 


(6^ — r^) + ba 


(32a) 

(32b) 


(33) 


and the exact solution in metric form is 


ds^ 


—dt^ + 





giAi exp(Ait)+g2A2 exp(A2t) 

hn/ 

(jj 

gi exp(Ait)+g 2 exp(A 2 q 



[dr'^ + r2 {dO"^ + sin2 9d(j)^')\ . 


(34) 

The new exact solution fl3T)l is similar, in structure, to the solution found by Goven- 
der and Thirukkanesh— . However our model results from a different physical scenario 
since the atmosphere of our star does not contain the cosmological constant, but com¬ 
prises a two-fluid system in which one of the components is a string fluid. This model 
corresponds to geodesic heat dissipation in a relativistic star when the string fluid in 
the stellar atmosphere is undergoing diffusion. The solution fl3T)) can now be used 
in the framework of irreversible causal and noncausal thermodynamics to study the 
temperature evolution of the radiating star. 
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3 Physical analysis 


In order to establish the role and influence of the null fluid energy density p at the 
surface S of the radiating star we probe the temporal behaviour of the fluid and 
gravitational variables in the generalised scenario with m = m{v,r) described in this 
instance by the exact solution fl34p . We compare this with the standard case with 
m = m{v) when p = 0, using the solution which can be found in Thirukkanesh and 
Maharaj— . We point out that these models are cast in natural units and consequently 
all physical parameters and mathematical constants are dimensionless. 

For the purpose of this investigation we plot the quantities B, p, z, m 

and ^ in Fig. 1-12. We show that the null density p drastically affects the physical 
behaviour of the model in the graphical plots. For the purpose of the graphs we make 
the following choice for the constants that appear in the exact solution: 

a = 1, 6 = 10, d = —1, gi = g2 = 1- 

For these values, flMll places a constraint on the allowed values for the energy density 
p and it becomes clear that 0 < p < 1. This restriction on p is unique for this model 
and our particular choice of values; it is therefore model dependant and will differ from 
another exact solution. 

3.1 Temporal evolution of the potential B 

The behaviour of the gravitational potentials i?p=o and Bp^Q is given in Fig. 1. It is 
evident that the potential B is always signihcantly larger in the case when p 7 ^ 0 than 
in the limiting scenario when p = 0. In the former, B is regular and decreases steadily 
with time, and increases in the latter. This trend suggests that the exterior null fluid 
has a suppressing effect on the gravitational held at the surface. We observe that this 
may in part be due to the nature and structure of our exact solution, and can possibly 
be somewhat different in another model. 
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Potential B 



Figure 1: Temporal evolution of the gravitational potential B at the stellar surface 

3.2 Temporal evolution of the fluid pressure py, 

Here we examine the temporal behavior of the radiating fluid pressure at the boundary 
of the object. The pressure, for the two cases, is dehned by the junction condition 
(IH]). Fig. 2 shows that for the case of pressureless radiation only, the fluid pressure ps, 
at the boundary is very small (p ps 4 x 10“^) and decreases monotonically with time 
until it is negligible. This is consistent with what has been theoretically established 
for stellar models in both the Newtonian and relativistic limits. On the other hand, in 
Figure 3 it can be seen that when p 7 ^ 0, p is slightly suppressed and becomes negative 
at the boundary. This is not surprising considering the general form of (0). Another 
reason that supports this result is that when p = 0 , the pressure is already close to 
zero, and if the density which is always strictly positive (however small in magnitude 
it may be) is taken into account it is clear that this will induce a further reduction. 
We have also investigated the effect of varying the magnitude of p and the family of 
prohles in Fig. 3 reveal that as the density is increased within the the allowed range, 
the pressure becomes increasingly more negative. This trend suggests that p scales 
inversely with p at the surface of the radiating star. It should be emphasized that 
despite this behaviour, the pressure is still only just slightly negative and does not in 
any way make the model physically unreasonable. 
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Figure 2: Temporal evolution of the radiating fluid pressure in the absence of p 
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Figure 3: Temporal evolution of the radiating fluid pressure in the presence of p with 
varying magnitudes 

3.3 Temporal evolution of the fluid luminosity Ly 

Utilising the dehnitions fll9ap and fllQbp along with the solution fl5T|l we are able to 
generate numerical prohles for the luminosity of the radiating fluid at the surface of 
the relativistic body and the luminosity as seen by a distant observer. In Fig. 4 we 
hnd that Ls evolves in an acceptable manner and decreases through positive values 
until it eventually becomes substantially low. The plot also indicates that the initial 
magnitude of the luminosity at the boundary is very small (Ls ~ 2.0 x 10“^), which may 
be justihable for stellar objects that are ultra-dense and possibly in the post neutron 
star phase of its evolution. After an extended but hnite period of time the star becomes 
virtually non-luminous and may not be observable. This is certainly consistent with 
the fact that the so called quark and quark-gluon stars have not yet been detected 
with current observing instruments. The family of prohles in Fig. 5 exhibits very 
interesting behaviour for Ls in the more general situation when p is non-vanishing. 
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Figure 4: Temporal evolution of the radiating fluid luminosity in the absence of p 
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Figure 5: Temporal evolution of the radiating fluid luminosity in the presence of p with 
varying magnitudes 

In general, the luminosity is suppressed and changes through negative values as the 
magnitude of p is increased; Ls at early times is less negative, and at late times tends 
to zero. The overall behaviour is similar to that found in Fig. 4 and again may be 
applicable to stars in the strong gravity regime. In Fig. 7, when there is only radiation 
in the exterior, the luminosity at inhnity is singular closer to the centre of the fluid 
distribution and evolves through positive but small values with r. It is also clear that 
for large r (r —)■ oo), L^o tends to zero. This suggests that a horizon will form at the 
end of the collapse of the radiating object. Fig. 6 on the other hand indicates that 
Loo, when p 7 ^ 0 , is hnite at r = 0 and evolves through negative values until eventually 
going to zero. Again, this behaviour demonstrates the formation of a horizon at late 
times. 
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Figure 6: Temporal profile of the radiating fluid luminosity at infinity in the presence 
of p 
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Figure 7: Temporal prohle of the radiating fluid luminosity at inhnity in the absence 
of p 

3.4 Temporal evolution of the gravitational redshift 

The plots in Fig. 8 show that Zp^Q > Zp=Q always. Furthermore, it is clear that the 
redshift, in both the standard Vaidya and generalised Vaidya models, evolves continu¬ 
ously and regularly through negative values. This is probably largely due to the nature 
of the exact solutions. 
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Figure 8 : Temporal evolution of the gravitational redshift at the stellar surface 

3.5 Mass flow across the boundary 

Here we investigate the temporal behaviour of the mass flow and rate of mass flow at 
the radiating surface. The numerical prohle in Fig. 9 illustrates that in the case of null 
dust (p = 0 ), the mass flow across the radiating surface decreases steadily with time, 
through positive values, until it terminates completely at later times. The situation 
when the flow of mass stops altogether corresponds to when the stellar object departs 
from its radiative phase and transits into an equilibrium state, usually towards the 
very late stages of gravitational collapse. For the model where the additional null 
fluid (p 7 ^ 0) is present in the exterior. Fig. 10 indicates that m increases steadily 
with time and eventually converges to a non-zero limiting value. At this point the 
star is still losing mass but at a uniform rate. Hence, we may conclude that in the 
standard Vaidya model the mass flow naturally comes to a stop. However, in the 
more realistic and generalised Vaidya case, mass continues to flow across the boundary 
for an extended period of time. This is strong motivation for the idea that the null 
fluid distribution in the local exterior is actually sourced from within the radiating 
body and is released along with radiation during dissipation. This must then mean 
that the null fluid is more likely to be a distribution of particles like neutrinos which 
is consistent with what is believed to be the case for most compact objects (see the 
earlier treatments of Misner— and Glass—). In Fig. 11 we observe that mass flow rate 
across the surface, decreases fairly rapidly and ultimately goes to zero when the flow 
ceases. The prohle in Fig. 12 indicates that for p 7 ^ 0 {dm/dt) is always negative and 
would seem to contradict the prohle in Fig. 10. This is probably due to the form of 
the gravitational potential in the solution fl3T)) and the choice of parameter values that 
have been considered. 
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Figure 9: Temporal evolution of the mass flow across the stellar surface, in the absence 
of p 



Figure 10: Temporal evolution of the mass flow across the stellar surface, in the presence 
of p 


3.6 Energy conditions with p ^ 0 

In any radiating relativistic model in astrophysics, a test of the strong, dominant 
and weak energy conditions are imperative in order to establish whether the model 
is valid and physically reasonable. This has been demonstrated in the seminal study 
by Kolassis et al— and in the recent treatment by Govender et al'^ , for dense stellar 
objects with heat conduction. The weak, dominant, and strong energy conditions are 
described respectively, in the generalised Vaidya spacetime, by the following 
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Figure 11: Temporal evolution of the mass flow rate across the stellar surface, in the 
absence of p 


(dm/dt)p^o 


0 


20 


40 

t 


60 


80 



Figure 12: Temporal evolution of the mass flow rate across the stellar surface, in the 
presence of p 


W = p — {qB — p) + A > 0, (35a) 

D = p — 3 {qB — p) + A > 0, (35b) 


S = {qB — p) + A > 0, (35c) 

where A is given by 

A = y/[(g5-p) + g]2_4g2. (36) 

In the limit when p = 0 the classical dehnitions for the energy conditions are regained. 

The corresponding numerical prohles for the parameters W, D, and S are given 
in Fig. 13-15. Fig. 13 demonstrates that the weak energy condition is satished since 
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0.050 


Figure 13: Temporal evolution of the weak energy condition parameter at the radiating 
surface, subject to the influence of the null fluid energy density p 



Figure 14: Temporal evolution of the dominant energy condition parameter at the 
radiating surface, subject to the influence of the null fluid energy density p 

the parameter W is regular, continuous and non-negative for all times. The dominant 
energy condition is also obeyed as seen from the prohle for D in Fig. 14. However, the 
strong energy condition is clearly violated since S is always non-positive, in Fig. 15. 
We must point out, though, that this result is not entirely detrimental for the model 
since it has been established by Kolassis et that the strong energy condition 
can be allowed to fail provided that the overall pressure of the matter distribution 
is negative. This works perfectly well with the generalised model with m = m{v,r) 
since the pressure P of the null fluid is negative for a constant density (as highlighted 
earlier in section 2.2) and in addition, the radiating fluid pressure is also negative (from 
Fig. 3), at the surface of the star. We may then conclude that the generalised null 
fluid model satishes collectively, all of the energy conditions and remains physically 
acceptable. 
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Figure 15: Temporal evolution of the strong energy condition parameter at the radiat¬ 
ing surface, subject to the influence of the null fluid energy density p 

3.7 The influence of the exterior null fluid on the nonadiabatic 
gravitational collapse of the relativistic star 

A relativistic star can radiate away its heat energy in two possible scenarios: in the 
process of nonadiabatic gravitational collapse and when the object is in hydrostatic 
equilibrium undergoing thermal cooling. Our exact solution can be applied to the for¬ 
mer case in particular and it would be a point of interest to determine what effect the 
exterior null fluid energy density p could have on the collapse rate at the boundary. 

The rate 0 of gravitational collapse for a spherically symmetrical shear-free fluid 
with geodesic particles is given by 


0 


35 

AB’ 


(37) 


for a comoving 4-velocity = 6q. With the solution (13T|) . for the situation where 

the additional null fluid is present on the outside, the analytical form for the rate of 
collapse is 


3 (6^ — r^) {b^(fp — a^) sech 


0p7^O - 


n 2 


¥\/ ^ ~ p 


2bd 


(V + r2) a + bd (V - r^) yj^ - p tanh \J - P 


(38) 


In the standard Vaidya model considered by Thirukkanesh and Maharaj— , with p 
we have that 


0 p=o — — 


3a (6^ — r^) exp j 


exp 


a(t+e) 

bd 




exp 


a(t+e) 

bd 


= 0 , 

(39) 
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Figure 16: The collapse rate for varying null fluid energy density 



Figure 17: The collapse rate for p = 0 

It is quite clear from equations fl38l) and fl39l) that in the general case 0p^o has a strong 
dependence on the density p and is far more involved than in the standard case 0p=o- 
The numerical prohles corresponding to fl38|l and (13911 are given respectively by Fig. 
16 and 17. 

From Fig. 16 it is evident that the rate of collapse is negative when the null fluid is 
present in the exterior of the radiating body. In other words the non-vanishing energy 
density p actually slows down the collapse process until such time when it reaches a 
complete halt and the relativistic star becomes stable again. It can also be seen that 
as the magnitude of the density is increased, 0p^o becomes less negative at earlier 
times and eventually goes to zero at late times. The plot in Fig. 17 illustrates that for 
vanishing density, the rate of collapse is always positive and indicates that the collapse 
itself is occurring reasonably faster. 
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4 Discussion 


In this paper we modeled a spherically symmetric relativistic radiating star that has 
vanishing shearing stresses and radial heat flow in the interior while its local exterior 
region consists of a two-flnid system comprising pressnreless radiation and a general 
nnll flnid. The Type II flnid possesses a non-zero pressnre and energy density and in 
particnlar we have investigated the inflnence of the density p on the physical behaviour 
of the dissipating stellar fluid at the surface. We demonstrated that the generalised 
junction condition p = qB — p for a null fluid scenario can be integrated as a nonlinear 
differential equation, at the boundary and generated an exact solution in the special 
case when the interior fluid particles are in geodesic motion and the exterior density 
is constant. Our solution, although similar in structure to that obtained by Govender 
and Thirukkanesh— , has a very different interpretation, depending strongly on the 
magnitude of p and can be cast in a more realistic astrophysical context. The gen¬ 
eralised two-fluid model with our exact solution has been compared to the standard 
model for pure radiation, studied by Thirukkanesh and Maharaj— , by constructing 
numerical prohles for various physical parameters associated with the stellar fluid at 
the boundary. 

We established that the constant null fluid energy density has a marked impact and 
signihcantly reduces the pressure ps and luminosity when compared to the limiting 
case when p = 0, where the pressure and luminosity are positive, hnite and consistent 
with well established results for extremely dense stellar objects in the strong gravity 
regime. The flow of mass across the radiating surface was also probed and it has been 
demonstrated that the density p actually enhances the outgoing flow; this is suggestive 
that the null fluid on the outside could more likely be interpreted as a distribution of 
particles that are sourced from nonadiabatic and particle production processes origi¬ 
nating within and at the surface of the star. This is in contrast to other treatments 
where the null fluid has been considered to be an ambient string system in four dimen¬ 
sions. In light of this our model may offer a more realistic description for dense stars in 
a relativistic setting. By utilising the interior matter variables and the exact solution 
we also tested the strong, dominant and weak energy conditions when the constant 
density null fluid is prevalent. Numerical prohles indicate that the weak and dominant 
conditions are satished while the strong condition is violated. However this does not 
severely affect the validity of the model in any way as the overall pressure is negative. 
Since radiating models are crucial for studying the nonadiabatic gravitational collapse 
of stellar objects, the rate of collapse, for both situations has also been examined. It 
has been found that p slows down the collapse process while in the case of null radiation 
only it occurs faster. Our model can be further improved on by including the effects 
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of shear and bulk viscosity and analysing the dynamical stability of the stellar matter 
conhguration. Furthermore, if the null fluid is to be taken as a system of particles 
like neutrinos for example then the model can be appropriately adapted by incorporat¬ 
ing non-vanishing neutrino fluxes in the interior. However, these considerations would 
consequently render the model more complicated and in particular the junction condi¬ 
tion increasingly difficult to solve, even with the robust mathematical techniques that 
are currently available. These endeavors will be the subject of an ongoing investigation. 
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